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We re-analyze the LHC bounds on light third generation squarks in Natural Supersymmetry, 
where the sparticles have masses inversely proportional to their leading-log contributions to the 
electroweak symmetry breaking scale. Higgsinos are the lightest supersymmetric particles; top and 
bottom squarks are the next-to-lightest sparticles that decay into both neutral and charged Higgsinos 
with well-defined branching ratios determined by Yukawa couplings and kinematics. The Higgsinos 
are nearly degenerate in mass, once the bino and wino masses are taken to their natural (heavy) 
values. We consider three scenarios for the stop and sbottom masses: (1) tii is light, (11) tL and 
bh are light, and (HI) tR, tL, and 6l are light. Dedicated stop searches are currently sensitive 
to Scenarios II and III, but not Scenario 1. Sbottom-motivated searches (26 -I- MET) impact both 
squark flavors due to i ^ bx^ as well as 6 — > fexi,2, constraining Scenarios I and HI with somewhat 
weaker constraints on Scenario II. The totality of these searches yield relatively strong constraints 
on Natural Supersymmetry. Two regions that remain are: (1) the "compressed wedge", where 
{mq — \fi\)/mq <^ 1, and (2) the "kinematic limit" region, where rUq ^ 600-750 GeV, at the kinematic 
limit of the LHC searches. We calculate the correlated predictions for Higgs physics, demonstrating 
that these regions lead to distinct predictions for the lightest Higgs couplings that are separable 
with ~ 10% measurements. We show that these conclusions remain largely unchanged once the 
MSSM is extended to the NMSSM in order to naturally obtain a large enough mass for the lightest 
Higgs boson consistent with LHC data. 



I. INTRODUCTION 

Natural Supersymmetry is the holy grail of beyond- 
the-standard model physics. It contains a sparticle spec- 
trum where sparticle masses are inversely proportional to 
their leading-log contributions to the electroweak symme- 
try breaking scale. At tree-level the electroweak sym- 
metry breaking scale is determined by balancing the 
Higgsino mass-squared against the scalar Higgs mass- 
squareds. This implies the leading contribution to elec- 
troweak symmetry breaking comes from the Higgsino 
mass itself, and thus implies the Higgsinos are the lightest 
sparticles in Natural Supersymmetry. The next largest 
contributions come from one-loop corrections from the 
stops. We consider the three scenarios: (I) in is light, 
(II) II and 6l are light, and (III) II, in, and &l are light. 
This spans the space of possibilities for various stop (and 
sbottom) mass hierarchies consistent with Natural Super- 
symmetry. After this comes the contributions come from 
the wino and gluino (in the MSSM), but their masses can 
be several times larger than the stop masses, given their 
comparatively suppressed contributions to electroweak 
symmetry breaking. Natural Supersymmetry suggests 
the lightest electroweakinos can be nearly pure Higgsino- 
like states. 

This spectrum is well-known plldOj. and the LHC ex- 
periments have already provided outstanding constraints 
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on simplified models involving light stops [¥T| - H5] . light 
sbottoms [m SS] , and gluinos that decay into these spar- 
ticles I46H49I . Further improvement in the bounds may 
be possible with specialized search strategies, for recent 
examples see [5DH5S]. However, the results presented thus 
far typically make strong assumptions about branching 
fractions [BR{ii ^ i^i) = 1 or BR{it ^ 6xf ) = llQ 
In addition, in cases where both a light chargino ana 
a light neutralino are present, the results assume cer- 
tain mass hierarchies: m-± =0.75mr -|-0.25m<;o 1431 or 
m~± — 2x m^o 141 or m-± —m^o > 50GeV HI]. These 
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assumptions make it difficult to extract the true bounds 
on Natural Supersymmetry. Consequently, we have un- 
dertaken a re-evaluation of the constraints on Natural 
Supersymmetry using the existing LHC results on sim- 
plified models involving light stops and sbottoms. 

It is also well-known that there is an intricate interplay 
between a light third generation and Higgs physics. Su- 
persymmetry predicts the mass of the lightest Higgs bo- 
son to high accuracy through radiative corrections that 
are dominated by just the third generation squarks [66l - 
[75] . If the third generation squarks are collectively light 



^ The notable exceptions are the two recent ATLAS searches for 
ii — > ix5 in the 1-lepton mode 41 and all-hadronic mode |42) . 
where constraints on the branching fraction BR{t\ — > tX\) were 
shown, assuming the remaining of the branching fraction is unob- 
servable. Natural Supersymmetry, however, predicts branching 
fractions into several channels that are observable, as we will see. 



(say, < 1 TeV), the predicted mass of the Higgs bo- 
son is too small to be compatible with the ATLAS and 
CMS observation [771 [75] of a 125 GeV Higgs-like bo- 
son (e.g. [Ml [371 Sni [TSHHI])- On the other hand, light 
third generation sparticles can significantly modify the 
detailed properties - production cross section and decay 
rates - of the lightest Higgs boson [5^105| . 

We consider the effects of Natural Supersymmetry on 
the detailed properties of the lightest Higgs boson. Here 
we are not interested in maximizing a particular decay 
channel or fitting to the existing Higgs results, but in- 
stead we endeavor to simply understand the characteris- 
tics that Natural Supersymmetry has on Higgs physics. 
Our main result is to overlay the modifications to the 
Higgs physics onto the allowed parameter space of Natu- 
ral Supersymmetry. Two interesting regions emerge. In 
the "compressed wedge" region where {mq — \fi\)/mq <C f 
and rriq can be small, the effects on Higgs physics are to 
enhance the inclusive (gluon-fusion dominated) cross sec- 
tion crJJJggjyj by 10-30% simultaneous with a slight reduc- 
tion of BRih — 77)mssm by up to 5%. By contrast, in 
the "kinematic limit" region where niq > 600-750 GeV, 
there is a slight enhancement of BR{h — > 77)mssm by up 
to 5%, with the inclusive (gluon-fusion dominated) cross 
section crJJJggj^ within a few % of the Standard Model re- 
sult. While the experimental situation the LHC collabo- 
rations is not yet settled, it is already clear that these two 
regions lead to distinctly different effects on Higgs prop- 
erties that can be probed with ~ 10% measurements. 

Given light stops and sbottoms, we must consider the 
supersymmetric prediction for the lightest Higgs boson 
mass. We assert that Natural Supersymmetry - in the 
MSSM - is simply incompatible with obtaining a lightest 
Higgs boson mass consistent with the LHC data. This 
point has been emphasized in some recent work, for ex- 
ample [8nil06lll07j . Hence, we do not restrict the third 
generation squark masses to obtain a given lightest Higgs 
boson mass. Instead, we assume there is another contri- 
bution to the quartic coupling that is sufficient to aug- 
ment the MSSM contributions, resulting in a Higgs mass 
that matches experiment, rrih ~ 125 GeV. Not specifying 
this contribution would seem to be fatal flaw of our anal- 
ysis. We show that simple extensions of the MSSM, in 
particular the next-to-minimal supersymmetric standard 
model (NMSSM), can give both a sufficient boost to the 
quartic coupling with negligible effects on the Higgsino 
mass spectrum and the decay chains that we consider 
here. Specific examples of NMSSM parameter choices 
that realize our assertion are given in Appendix \K\ 

We do not consider the gluino in this paper. The 
gluino contributions to the electroweak symmetry break- 
ing scale may be significant in the MSSM, given the 
existing searches that suggest the gluino must be heavier 
than 1-1.3 TeV, depending on the search strategy [i51 - 
l49l I108j . However, the size of the gluino contribution to 
electroweak symmetry breaking is model-dependent: A 
Dirac gluino has a substantially smaller contribution to 
the electroweak symmetry breaking scale compared with 



a Majorana gluino, when the leading-log enhancements 
are included, allowing a Dirac gluino to be substantially 
heavier |109H111] . In addition, the search strategies for 
a gluino depend on its Majorana or Dirac character. 
One of the most important search strategies - involving 
same-sign dileptons (such as [46l|49]) does not provide a 
constraint on a Dirac gluino. 



II. MASS HIERARCHY IN NATURAL 
SUPERSYMMETRY 

A. Contributions to the EIectroweal< Scale 

In the minimal supersymmetric standard model 
(MSSM) the electroweak symmetry breaking scale is de- 
termined by, at tree-level |112j . 



tan2/3 + lm^ -m^ 1 1 

^^^^FJ^l 2 2"^«„-2-H.-H . 

(1) 

In saying "contribution to the electroweak scale", it is 
understood that the supersymmetric and supersymme- 
try breaking parameters are adjusted to obtain the value 
already determined by experiment. Here we are inter- 
ested in the relative size of and the loop corrections 
to the electroweak breaking scale, i.e., Mz- 

For tan/3 very near 1, the coefficient of the first term in 
Eq. ([T]) becomes large, because the _D-flat direction in the 
scalar potential is not lifted, and thus implies increased 
sensitivity to the supersymmetric parameters. The sen- 
sitivity is most easily understood by eliminating depen- 
dence on rn?u^ using the tree-level relation [112] 



'Ha 



to obtain 



-Mi 



1 



tan^ 13 - 1 



tan^ /3 + 1 
tan^ /3 



At large tan/3, however, Eq. (|3| simplifies to 
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and eliminates dependence on m\. Generally, we have 
taken tan/? ~ 10 for the analyses to follow, and thus 
the heavy Higgs scalars that acquire masses near niA can 
be readily decoupled from our analysis. However, the 
smaller tan/3 region reappears in our discussion of the 
NMSSM in Appendix |^ where the the relative contribu- 
tions to the electroweak symmetry breaking scale become 
more complicated for the NMSSM scalar potential. 

With Eq. Q in mind, we can compare the relative 
importance of different contributions to the electroweak 
symmetry breaking scale by normalizing to M'^/2 [113] 
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The tree-level and one-loop contributions are well-known (e.g., [371 lll^l ) 
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In the MSSM, there are also important one-loop contributions from a Majorana wino and two-loop contributions from 
a Majorana gluino (e.g., |57llll2) ) 
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Here we somewhat arbitrarily chose to normalize all of 
our numerical evaluations to a factor of 10 times Af|/2, 
as well as normalizing the size of the leading-logs to 
Amoss/w = 20|^ This small ratio implicitly assumes a 
low scale for the messenger sector, and thus the small- 
est sensitivity of supersymmetry breaking parameters to 
the electroweak breaking scale. This provides suggestive 
values for \fj,\, the stop masses, and in the MSSM, the 
wino and gluino masses. As these parameters signifi- 
cantly differ from these suggestive values, their relative 
importance to determining (or fine-tuning to determine) 
the electroweak scale is altered accordingly. In particular, 
we see that = 200 GeV gives a comparable contribu- 
tion to a pair of stops at — — 450 GeV. 

The Natural Supersymmetry predictions for the wino 
and gluino mass depend on whether they acquire Majo- 
rana or Dirac masses. Already we see that if the wino 
acquires a Majorana mass, its mass is expected to be 
nearly 1 TeV. A Dirac wino would have a mass consider- 
ably larger. Similarly a Majorana gluino is expected to 
be 1.2 TeV (with the normalization of the logs as given 
above), and again significantly larger than this if it ac- 
quires a Dirac mass. 

For the purposes of this paper, we assume the gluino is 
either sufficiently heavy so as to not lead to collider con- 
straints (in practice, this means a Majorana gluino needs 



Except for Anicss/IAfsl = \/20, since a conservative interpreta- 
tion of LHC bounds is tliat the gluino already exceeds 1.3 TeV 
in viable scenarios. 



to be above about 1.3 TeV [46ti49l fT08] ) . or it acquires 
a Dirac mass, in which case its Natural Supersymmetry 
mass is well out of range of the LHC. We assume the 
wino and bino acquires ~ 1 TeV masses, but our results 
are largely insensitive to this choice. 



B. Higgsino mass splitting 



In the limit Mi ^2 ^ ImIj'^? the lightest chargino and 
the lightest two neutralinos are Higgsino-like and nearly 
degenerate in mass. The leading contributions to the 
mass difference at order 1/Mi_2, 
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FIG. 1. The three scenarios of Natural Supersymmetry considered in this paper. Scenario 1,11, and III are illustrated in 
the left, middle, and right panels. In our study we varied /i ~ m^o between 100 — 500 GeV and lightest stop masses mi_^ ~ 
250 - 1000 GeV. 



which we can write as 
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Clearly, the mass differences among the Higgsinos are 
just a few GeV when Mi^2 take on natural (heavy) val- 
ues. The decays xJjX^ ~^ Xi thus yield unobservably 
small energy in the decay products. The mass difference 
is, however, large enough that the decay rates are prompt 
on collider time scales and thus there are no macroscopic 
signatures in the detector (at least for wino and bino 
masses that do not far exceed ~ 1 TeV). Hence, for the 
purposes of LHC detection, , x\ 2 behave as neutral 
lightest supersymmetric particles that escape the detec- 
tor as missing energy. 



C. Simplified Models of Natural Supersymmetry 

Evidently from Eq. ([t]), the Natural Supersymme- 
try prediction for the stop masses depends on the sum 
TOq^ + '^^3 + l^tP- All other things considered equal. 
At ^ Q implies the sum juq + tti'^,^ must be correspond- 



ingly smaller to hold the sum + + \At\ con- 
stant. We therefore take At to vanish. While this might 
give some readers pause, regarding the stop contribu- 
tions to the lightest Higgs mass, recall that we have al- 
ready asserted that the MSSM is incapable of providing 
a sufficient contribution, and so the choice ^4^ = is not 
inconsistent with our approach. Instead, we consider the 



following mass hierarchies ("Scenarios") for Natural Su- 
persymmetry: (I) tfl is light, (II) tL and 6l are light, and 
(III) in, II, and are light. These scenarios span the 
space of possibilities for the stop (and sbottom) masses. 
We illustrate these scenarios in Fig.[l] The resulting mass 
eigenstates are given by |112j 

Scenario I ^\ — "^L + ™t + ^ur (12) 
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where = (T, - sin^ Sh') cos(2/3)M|. In Sce- 
nario III, we take the soft masses to be equal JTig^ — 
mcigj^ Since cos(2/3) < for tan/3 > 1, this im- 
plies Aufi > whereas A^^ < 0, causing ii ~ 
and t2 — Ir. Given that we specify soft masses, 
hi is always lighter than ii in Scenarios II and III. 
The mass difference is (50, 30, 20) GeV for m^^ — 
(200,400,600) GeV, corresponding to [mi^.mi) ~ 
[(260, 210), (435, 405), (620, 600)] GeV. Finally, we also 
impose — m^o > 50 GeV for reasons related to 

the details of the search strategies employed by ATLAS 
and CMS. We discuss this in the next Section. 

All other gauginos, all sleptons, and the first and 
second generation squarks are taken to be sufficiently 
heavy that they do not play a role in the low energy 
phenomenology, consistent with Natural Supersymmetry. 
We emphasize that the difference between Scenarios I, II, 
and III is not that the stops are believed to be far differ- 
ent in mass, but simply different enough in mass that the 



^ We also include the left-right squark mixing contribution 
mt/i/ tan /9, but this is suppressed by our choice of tan /3 = 10 as 
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low energy phenomenology is dominated by one of these 
three scenarios. 



III. COLLIDER BOUNDS ON NATURAL 
SUPERSYMMETRY 

A. Collider study setup 

The inputs for the spectra are the soft masses 
TOQ3,m„3, 771^3, the /i term, tan/?, the bottom and top 
Yukawas, and the weak scale v. For simplicity, we will 
assume the At, A}, terms are zero. However, by inter- 
polating between our results for the three given Sce- 
narios, it is possible to reconstruct qualitatively what 
happens when At^b 7^ 0. We set tan /3 = 10, and any 
"decoupled" particle in a given Scenario is, for purely 
practical reasons, taken to have mass 5TeV. Our re- 
gion of interest is m^o > 100 GeV, mj^ > 250 GeV. 
The lower bound on the LSP mass comes from the LEP 
bound on charginos ~ we must obey this bound since 
m^o ^ ^ '"^x^ ~ Finally, while the viability of stops 
with mass < 250 GeV remains an interesting question, 
we concentrate on mj^ > 250 GeV, consistent with the 
Natural Supersymmetry spectrum. 

We also impose an additional restriction on the pa- 
rameter space, namely, to not let the mass difference be- 
tween the squarks and the Higgsino become too small. 
The experimental analyses on stop production and decay 
through t — >• txi restricted m^^ ~"^x? > l'''^ GeV for AT- 



LAS semi-leptonic and all-hadronic searches [4TJ |42] , and 
> 200 GeV for the CMS semi-leptonic search gS]. Stop 
decays in compressed spectra often lead to to multiple- 
body final states that are difficult to model without bet- 
ter tools. Additionally, the collider limits on nearly 
degenerate spectra become sensitive to how the addi- 
tional radiation in the event (ISR) is modeled. We chose 
to simulate the sensitivity of these searches in Natural 
Supersymmetry for somewhat smaller mass differences, 
m/i - — m^o > 50 GeV. As we will see, we find the 
existing LHC searches are sensitive to Natural Super- 
symmetry with splittings this small. However this region 
needs to be interpreted with some care, since we are ob- 
taining constraints from both stop searches and sbottoms 
searches. Sbottom searches have somewhat different re- 
strictions on the kinematics, but since we chose a mini- 
mum mass difference between the soft mass toa and the 
lightest neutralino, the highly compressed region with 
respect to the sbottom and neutralino is not simulated. 
We therefore do not anticipate significant changes in the 
bounds for the parameter space we consider. 

For each scenario at a given mass point (p, "m-q), 
events are generated using PYTHIA6.4 |114j . We use 
CTEQ6L1 [115] parton distribution functions and take 
all underlying event and multiple interaction parameters 
to their values specified in Rcf . |108( I116j . The cross sec- 
tion is calculated by summing the next-to-leading order 
plus next-to-leading log (NLO + NLL) values [TTTHT^ 



using jl24j over all light (< TeV) 3rd generation sparti- 
cles in the spectrum. Following PYTHIA generation, the 
events are fed into DELPHES |125| to incorporate detec- 
tor geometry and resolution effects. We use the default 
DELPHES ATLAS and CMS detector descriptions, but 
modify the jet definitions to agree with the correspond- 
ing experiment: anti-fc^ algorithm, with size R = 0.4 for 
ATLAS and i? = 0.5 for CMS analyses. Additionally, 
while the experimental flavor-tag/fake-rates slightly dif- 
fer from analysis to analysis, we used a fixed 70% tag 
rate for all 6-jets that lie in the tracker {\r]j\ < 2.5). 

The signal simulation is used to derive the efficiency 
- the survival rate in a given bin of a particular anal- 
ysis. Given our simplified supersymmetry spectra, this 
efficiency is a function of the squark and LSP masses 
alone, i.e. for bin i we find ei{'m^_^,m^o). The product of 
the derived efficiency with the luminosity and the cross 
section (at NLO + NLL) is the number of signal events. 
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We derive exclusion limits by comparing s,; at a given 
{rhi,mj^o) with the number of signal events allowed at 
95% CL calculated with a likelihood-ratio test statistic. 
Specifically, the 95% CL limit on the number of signal 
events, Si 95 is the solution to 



0.05 



H,iPois(nj|6j -I- 5^,95) 
UiPois{ni\bi) 



(15) 



where the Ui is the number of observed events in a chan- 
nel «, bi is the number of expected SM background events, 
and we take the product over all orthogonal channels. We 
take both and hi directly from the experimental pa- 
pers. 

To incorporate systematic uncertainties, the number of 
background events in a bin is allowed to fluctuate: bi — 
bi (1 -|- Sbi). After multiplying by a Gaussian weighting 
factor, we integrate over Sbi, following Ref. |126j . We 
take the width of the Gaussian weighting factor to be the 
relative systematic uncertainty in a given bin quoted by 
the experiment jfj*, using the larger error if asymmetric 
errors are giverrl 

Pois{ni\bi + 95) -> 

6bi Gaus(J6j, )Pois(ni|6j (1 -|- 6bi) + 3^,95) 

(16) 

One may ask if a likelihood-ratio analysis is really 
needed, instead of just a rescaling of existing bounds. If 
the signal yield according to ATLAS/CMS was given for 
each (mj^ , m^o ) bin, then we could rescale and determine 



* The Gaussian integration is truncated such that the number of 
background events is always positive. 
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the yield, and thereby the exclusion bounds, in each of 
our Scenarios. However, such detailed information is not 
public and only the signal yields at specific benchmark 
points are given. In order to extrapolate yields away from 
the benchmarks, some model is needed, and for that we 
rely on the simulation method described above. 

Before describing the details of the searches we con- 
sider, it is important to emphasize that the absolute 
bounds we present are only approximate. To derive 
the signal efficiency we have used fast-simulation tools 
(DELPHES) whose energy smearing and tagging func- 
tions are approximations - usually optimistic - of the full 
detector effects. In multi-jet, especially multi-6-jet final 
states, the diff'erences between the fast and full-detector 
simulations add up, making it tricky for us to match the 
quoted absolute bounds on a given scenario. To improve 
the accuracy of the absolute bounds, the Scenarios pre- 
sented here should be studied by CMS/ATLAS them- 
selves, either as a dedicated reanalysis or using a tool 
such as RECAST |127| . Meanwhile, the relative bounds, 
i.e. the difference between Scenario I and Scenario II, are 
robust. 



B. Direct stop searches 

In this section we present the limits on the super- 
natural scenarios from the most recent LHC direct stop 
searches pP-'33' . Bounds from these searches are usually 
(though not always) cast in term of stops that decay ei- 
ther 100% of the time to a top quark and a neutralino 
or 100% of the time to a bottom quark and a chargino. 
Stops are searched for in several different final states, and 
the first two stop analyses we consider are semileptonic 
searches. While the details differ between the ATLAS 
and CMS searches (see Appendix [B| for the full analyses 
description), both require a hard lepton, significant miss- 
ing energy, and at least four jets, one of which must be 
tagged as a 6-jet. 

Running our three Scenarios through the CMS direct 
stop search [43j . we find the following exclusion contours 
in the (m^^ , m^o ) plane (Fig. j2j) . This is somewhat an 
abuse of notation - the horizontal axis actually corre- 
sponds to the mass of the lightest stop eigenstate for a 
given spectrej^ For comparison, we include limits from 
two "default" spectra (calculated in the same manner as 
our three Scenarios): 

1. Stop production and decay with 100% branching 
fraction to top plus neutralino. The decay is carried 
out using phase space alone, so the decay products 
are completely unpolarized. This setup is exactly 
the CMS simplified model T2tt [108]. 



2. Right-handed stop production followed by decay to 
a bino-like LSP plus a top quark. In practice we 
take the exact setup for Scenario I but replace swap 
the roles of /i and Mi . This spectrum is close to the 
default signal model used by ATLAS. As the hand- 
edness of the stop and the identity of the LSP are 
fixed, the polarization of the emerging top quark is 
also fixed. 

By comparing our Scenarios with the stop signal models 
usually used, we can see how the Higgsino-like nature of 
the LSP and the hierarchy of the third generation squarks 
effects what regions of parameter space are allowed. The 
comparisons also give some indication of how well our 
simple analysis matches the full ATLAS/CMS results. 

We can understand the strength of the bounds on our 
Scenarios by looking at the branching ratios and final 
states of our spectra. As we have decoupled the gaug- 
inos in all of our setups, the decays of the stops and 
sbottoms are governed entirely by the Yukawa couplings. 
For example, in Scenario I, all decays come from the 
top- Yukawa; up to kinematics, this yields a 50-50 split 
between decays to top quark plus neutralino and bot- 
tom quark plus chargincj^ Due to the degeneracy of 
the chargino-neutralino sector, chargino decay products 
are all extremely soft. In particular, the leptons from a 
chargino decay are far too soft to trigger the analysis re- 
quirements for the stop analysis, thus the only source 
of leptons is from the stops that decay to a leptoni- 
cally decaying top quark. Additionally, mixed decays 
iiil t{-^ ivb) + + b + Xi may have a hard lep- 
ton, but they typically have fewer jets than required for 
a stop analysis. Therefore, only the fraction of events 
where both the stop and antistop decay to top -I- neu- 
tralino have a high probability of passing the analysis 
requirements. As an final suppression, because the light 
stop in Scenario I is (almost) entirely right handed, the 
top quarks it yields are left-handed. Due to the V-A 
nature of the weak interaction, left-handed tops have a 
softer lepton spectrum, and thus the leptons that the 
stops decays do create are less likely to pass the analysis 
cuts |128) . The combined effect of these suppression fac- 
tors is that there is no bound from the CMS direct stop 
search on Scenario I. 

Similar logic works to understand the bounds on Sce- 
nario II and III. In Scenario II, both the 6^ and are 
produced. Up to effects of 0(j/f, tan /3) and ignoring kine- 
matics, the only decay channel possible is t ^ t + for 
the stop and 6 — > t 4- fo^' the sbottom. The stop 
therefore decays in exactly the same way as in the de- 
fault scenario, so we expect the bound to be at least as 
strong as the T2tt bound (with the added effect that the 



For example, in Scenario III the spectrum also contains a second 
stop and a sbottom - all three states are produced and analyzed 
when deriving the analysis efficiencies, though limits are still 
placed in terms of the lightest stop eigenstate 



The neutralino branching fraction is further split: roughly 50% 
to X\ S'^d 50% to However, this distinction does not make 
affect our analysis, since the two states have essentially the same 
mass 
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FIG. 2. Limits on the various 3rd generation scenarios com- 
ing from the CMS direct stop search (semileptonic channel). 
The a;-axis corresponds to the lightest physical stop mass t\ 
in each Scenario. In Scenarios II and 111, one 6-squark is 
present with a physical mass that is slightly lighter, given by 
Eq. ( 13 1. The orange contour shows the 95% exclusion bound 



on Scenario III, the green contour is the bound on Scenario 
II, and there is no bound on Scenario 1. The black, dashed 
contour is the bound from this analysis using CMS simplified 
model T2tt that involves direct production of stops which de- 
cay solely to unpolarized top plus neutralino; ti — ^ t + Xi- 
The brown dashed line shows the limit on a second default 
scenario: right-handed stops decaying to top plus bino. The 
difference between the black and brown dashed lines gives 
some indication how the polarization of the top can affect 
limits. The black dotted line is mj — m^o = 150 GeV, which 
is the self-imposed restriction on the CMS analysis, since ISR 
is not properly taken into account when the signal is gener- 
ated with PYTHIA. We have also restricted our re-analysis in 
the "compressed wedge" region [where {niq — m^o ) /rriq <gi 1] , 
requiring ?tIq^ — m^o > 50 GeV, that results in the excluded 
region extending slightly to the left of the black dotted line. 
See the text for details. 



top is always right-handed and thus the emitted lepton is 
harder than in the unpolarized case). The bound is actu- 
ally stronger because the sbottom decays also contribute; 
the chargino in a sbottom decay is indistinguishable from 
a neutralino, so the final state from a sbottom decay is 
virtually identical to the stop case. The only place the 
bound on Scenario II may weaken is close to or below 
the i + x° threshold, where t ^ b -\- xt decays become 
important. Finally, we expect an even stronger bound in 
Scenario III. In addition to the h decay that contributes 
exactly as in Scenario III, there are now two stop states 
and both states will contribute to the stop search. These 
suspicions are confirmed in Fig. [2j 



FIG. 3. Limits on the various 3rd generation scenarios coming 
from the ATLAS direct stop search (semileptonic channel). 
Contours are the same as in Fig. [2] 



Moving to the ATLAS semi-leptonic stop search, we 
find similar results, shown in Fig. |3] This is not surpris- 
ing as the search criteria are very similar to the CMS 
stop search - a single hard lepton and four or more 
hard jets. The biggest difference between the ATLAS 
and CMS semi-leptonic stop searches is that ATLAS re- 
quires a "hadronic top candidate" - a three-jet subsys- 
tem with mass between 130 GeV and 205 GeV - in all 
events. This requirement, along with slight changes in 
the analysis variables (see Appendix |b| lead to differ- 
ent limits, but the qualitative message is the same as in 
the previous case: scenarios with ^ <C , m^;^ are 
not bounded by these searches because the stops decay 
preferentially to 6 + x J and therefore lack sufficient hard 
leptons and jet multiplicity, while scenarios with light 
ti, 6l are bounded tighter than the benchmark t ^ t+}(^ 
scenario because both the stop and the sbottom decays 
contain top quark^ 



Comparing our bounds for t\ — > tx*^ (T2tt model) with the 
exclusions from ATLAS, we see a discrepancy - our bounds 
are weaker by O(IOO) GeV. The fact that the DELPHES-based 
bound is different from the quoted number is not surprising, but 
the discrepancy is somewhat larger than expected. ATLAS has 
provided a cut-flow, at least for some benchmark (mfj^i^^o) 
points, which allows us to pinpoint the difference to the rrijjj 
cut (relative efficiencies of cuts either before or after this cut 
match to 0(10%)). We suspect the reason the rrijjj cut is dis- 
crepant is that the jet-energy resolution in DELPHES is overly 
optimistic. If the jets retain too much of their energy, then the 
rajjj distribution will be shifted to higher values (relative to the 
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bounds we find are shown below in Fig. [4j The trend of 
these bounds is similar to what we found in the previous 
stop searches, though the reasoning is slightly different. 
The bounds on Scenario II and III are nearly identical 
and rule out stops below 750 GeV for /i — 100 GeV. 
There is no significant bound on Scenario I due to the 
high fraction of decays to b + Xi] stop decays to bot- 
tom quarks do not contain enough hadronic activity to 
efficiently pass the jet multiplicity cuts in this analysis. 

Summarizing the direct stop searches, scenarios with 
degenerate, light Higgsinos and Ir <C ^Lj^l are very 
weakly bounded, while the bounds on scenarios with light 
iL,6i, are quite tight, typically 100 GeV stronger than 
the bounds on the benchmark ix t + setup. Be- 
cause the direct stop searches are so insensitive to light 
tfi (with light /i), the bounds on Scenario II {1^,6^ and 
tfi all light) and Scenario III (only tL,bL light) are nearly 
identical. However, before we can draw any firm conclu- 
sions on Natural Supersymmetry, we must also consider 
the CMS and ATLAS searches tailored towards the de- 
tection of sbottoms. 



FIG. 4. Limits on the various 3rd generation scenarios coming 
from the ATLAS direct stop search (all-hadronic channel). 
Contours are the same as in Fig. [2] 

The final direct-stop analysis we explore is an all- 
hadronic search performed by ATLAS using 20.5 fb"^ of 
datgj^ Unlike the previous stop analysis, no leptons are 
involved. Instead, stops are searched for in events with 
multiple hard jets (6 or more), at least two 6-jets, and 
substantial missing energy. To suppress multi-jet QCD 
backgrounds, the jets in the event are required to form 
two top-candidates - three-jet subsystems with invariant 
mass between 80 GeV — 270 GeV. When interpreted in 
terms of the ti — >■ t -I- x° benchmark scenario, ATLAS 
finds the strongest stop bound to date, nearly 700 GeV 
for massless neutralino. 

Applying these analyses to our three Scenarios, the 



full detector) and lost once the cut rrijjj < 205 GeV is imposed. 
If we increase the upper rrijjj cut by ~ 50 GeV, the signal ef- 
ficiency at the benchmark point agrees better with the quoted 
value, however this artificial shift will have uncontrollable impli- 
cations in the rest of the (m^^ , m-o) efficiency plane. Therefore, 
we stay with the quoted cuts and emphasize that the relative 
bounds between models are the most relevant. The strong sen- 
sitivity of the bounds to rrijjj also serves as a warning to the 
experiments since rrijjj is susceptible to effects from ISR, the 
underlying event, and pileup. 

ATLAS has performed a stop search in the dilepton final 
state 11291 assuming _BR[ti — > fex^] ~ 100% and using a va- 
riety of chargino-neutralino mass splittings (though none consis- 
tent with fj, <JC Ml, M2). As the search requires two leptons, the 
same issues raised for Scenario I will be present and we expect 
no bound. For Scenarios II and III, we expect stronger limits 
from the semi-leptonic search since the decay ii — > is dom- 
inant. For these reasons we do not explore the limits from the 
dileptonic searches on Natural Supersymmetry. 



C. Direct sbottom searches 

In Natural Supersymmetry, the stops can decay into 
b + xj*^, and thus dedicated searches for 6-quarks plus 
missing energy are vital to our analysis. In addition, in 
both Scenarios II and III, &i is present in the spectrum 
with rrij,^ ~ m^^ determined by uiq^. In this section 
we use the ATLAS and CMS searches that target direct 
sbottoms [in [45], since these studies focus on 6-jets and 
missing energy and are therefore independent of the mass 
splittings in the chargino/neutralino sector. 

To isolate signal-rich regions from background, AT- 
LAS/CMS sbottom searches require multiple high-px jets 
along with one or more flavor tags. Events with leptons 
are vetoed as a way to remove some tt background (the 
leptonic events are retained as control samples). More 
elaborate cuts are applied to further enhance the sig- 
nal depending on the collaboration and the target sig- 
nal mass. The default signal we will compare to is pair- 
production of sbottoms that decay solely to b quarks plus 
neutralino, a, bb + fx final state. Since it is identical to 
the CMS default signal model, we will refer to the default 
as T2bb as they do. 

To bound bb + signals, CMS [45] retains events 
with 2-3 jets and 1 or 2 5-tags. The visible objects in the 
event are partitioned into two "mega-jets". The degree 
to which these mega-jets balance each other, described 
with the ax variable [1301 1131j , as well as the net Hj- are 
used to further isolate the signal from background. The 
bounds from this analysis on the CMS default model and 
on our three scenarios are shown below in Fig. [5| 

The first thing to notice is that the sbottom search 
places a strong bound on Scenario I - roughly 'm^^ > 
600 GeV for m^o ^ ^ ^ 100 GeV and decreasing slightly 

as m^o is raised. 

^1 
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FIG. 5. Limits on our 3rd generation scenarios from the CMS 
direct sbottom search. Contours are the same as in Fig. [2] 
however there is now a bound, indicated in blue, on Scenario I. 
We have also added the bound (red dashed line) derived from 
applying this analysis to the T2bb simplified model, direct 
production of right-handed sbottoms with 100% branching 
fraction to a bottom quark and a neutralino. The remaining 
contours are the same as in Fig. [2] 

The bounds on Scenario I are weaker than the bounds 
on the T2bb scenario. This is because Scenario I yields 
more leptons - coming, as before, from stop decays to lep- 
tonic tops - so events from Scenario I are more likely to be 
vetoed. Also, the average number of jets is higher, push- 
ing the signal into jet bins not considered in the sbottom 
analysis. The same two effects also explain the difference 
in bounds between Scenarios I and II. In Scenario II, pro- 
vided TOj^ ^ mt -\- m^o , both stop and sbottom decays 
result in top quarks. The only di-top quark events that 
cleanly mock the bb + fx signal are fully leptonic events 
where both leptons are lost or lie outside the tracking 
volume. In all other events there is either a lepton or a 
larger jet multiplicity and the event is either vetoed or 
populates a region not usually considered as signal. The 
caveat to this argument is when m^^ < nit + "T-j^o . In this 

region, kinematics suppresses the ii — t- < + x'' mode and 
the (otherwise Yukawa suppressed) ii ^ b + mode 
becomes important. Decays to 6 -I- Xi^ are efficiently se- 
lected by the CMS search, explaining why the bound on 
Scenario II gets stronger the closer the stop mass gets 
to rrit + TO^o . The bound in the threshold region of Sce- 
nario II is actually stronger than in Scenario I since both 
t L and bL are produced and both decay to 6 -|- x when 
m^^ ~ ^mt+ m^o . As expected, the bound on Sce- 
nario III is the strongest and resembles the sum of the 
bounds on Scenario I and II. 




400 600 800 1000 

nij^ (GeV) 

FIG. 6. Limits on our 3rd generation scenarios from the AT- 
LAS direct sbottom search. Contours are the same as in 
Fig. [5] 



The ATLAS direct sbottom [44] search targets the 
same final state, bb + fx as the CMS search. However, 
the ATLAS search is more optimized to the topology with 
exactly two bottom jets, missing energy, and little other 
hadronic activity. A third hard jet is vetoed in the ma- 
jority of the analysis channels, and no channel tolerates 
4 or more jets. As a result, the ATLAS sbottom search 
is less flexible and not as well suited to events that con- 
tain top quarks. The bounds from the ATLAS sbottom 
search cast in term of our scenarios and the benchmark 
T2bb model are shown below in Fig. ?'iF] 



D. Combined Bounds 

Combining the three stop searches and two sbot- 
tom searches by taking the strongest limit at a given 
(r7ij^,m^o) point, we get the net excluded region for the 
three Scenarios. The excluded regions are displayed be- 
low in Fig. [T] along with the analogous regions for the 
default spectra. 



In the ATLAS analysis 11321 the sbottom search technique was 
used to constrain stop production, exactly as we are advocating 
here. However in that analysis, BR(tt — bxi) = 1 was assumed, 
so bounds presented there do not constrain scenarios with a light 
Higgsino, a key ingredient in Natural Supersymmetry. 
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FIG. 7. Limits on our 3rd generation scenarios from com- 
bining all CMS and ATLAS sbottom/stop searches search. 
Contours are the same as in Fig. [5] 



IV. IMPLICATIONS FOR THE HIGGS SECTOR 

In this section we study the implications of Scenario 
I, II and III on the supersymmetric Higgs sector. In 
supersymmetry the additional charged and colored de- 
grees of freedom can significantly modify the production 
cross section and branching ratios of the lightest (stan- 
dard model-like) Higgs boson [92Ul05j . Given the recent 
discovery of a particle consistent with a Higgs boson at 
rrih — 125 GeV [771 ESj) the modifications due to the 
additional charged and colored degrees of freedom have 
been extensively studied (40l HUHl [Ml HH [l33] . 



First let us consider the Higgs boson branching ratios. 
When the stop contributions are included, the modifica- 
tion to the decay rate of the Higgs boson into gluons is 
given by [l02l ITM] 



pMSSM 

ggh 

pSM 



1 



A/2iTt) ^ 
2 



1 



E 



(17) 



where v4i/2(^o) are the standard fermion (scalar) loop 
functions (e.g. [134] ). and = m\/4mf. Here mj^ > 
'nij_^, Of is the sfermion the mixing angle, and the cou- 
plings are given by 



"^1 



m^- / /If 2 ' 



fi 



2 I ^ h (18) 



in the decoupling limit. Hence in the limit of small 
mixing. Of ^ 0, the squarks enhance the decay rate of 
the Higgs boson into gluons. Similarly, for sbottoms 
the contributions are typically small except in the large 
tan /3 regime where sbottom contributions will interfere 
destructively with the top contribution. 

Light stops, sbottoms, staus and charginos also affect 
the decay of the Higgs boson into photons. Enhancing 
the decay rate of the Higgs to gluons due to light stops 
will lead to a suppressed decay rate of the Higgs into pho- 
tons due to the stop contribution destructively interfering 
with the ly-boson contribution (the dominant standard 
model contribution), while a light sbottom has the op- 
posite effect. Furthermore, depending on the sign of ft, 
a light Higgsino close to the LEP bound [135] can either 
enhance or suppress the photon rate |136j . Expanding 
in terms of the I/M2, where M2 is the Wino mass pa- 
rameter, we find the amplitude of the lightest chargino 
is 



2M^ 

\M2\m^o 



|Cq+^|Ai/2(t^o) 



(19) 



where Mw is the W-boson mass, M2 is the Wino mass, 
m^o ~ ij, and a is the mixing angle of the CP-even Higgs 
bosons. In the decoupling limit Ca+p ~ S2/3, this contri- 
bution becomes suppressed for large tan /3 |136j . There- 
fore the charged contributions are also relatively sup- 
pressed in an MSSM-like framework. 

In Fig.[8|we show the impact of the spectra in Scenario 
I, II and HI on modifications to BR{h — > 77), ainci, and 
ainci X BR[h — >■ 77). This is the principle result of our 
paper. We have taken /i > and — 00, however 
the results are nearly identical (to within ~ 5%) when 
Ma — 1000 GeV{^ Here we have assumed there is an 
additional contribution to the quartic coupling, raising 
the Higgs mass up to the experimentally measured value 
ruh ^ 125 GeV, such as the NMSSM-hke scenario de- 
scribed in Appendix [A[ 

We see that the modifications to the inclusive produc- 
tion cross section (dominated by gluon fusion) are at the 
10-30% level in the "compressed wedge", while rather 
small < 5% at the "kinematic limit" . The BR(h — 77) 
receives considerably smaller effects, between —5% to 
-1-5% across the parameter space of interest. These devi- 
ations are not large enough to be directly constrained 
by the measurements from the LHC [137] and Teva- 
tron |138j . however as we measure the Higgs production 
and branching ratios more precisely, we expect these de- 
viations to be observable at the LHC. 



For lower values of m^, the increased mixing between the two 
CP-even Higgs bosons leads to a slight further suppression in the 
branching ratios. 
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FIG. 8. Modifications to Higgs production and branching ratios in the decouphng hmit where Ma oo,At = Ab = G and 
M2 = 1 TeV. We have overlayed the direct search constraints found in Fig. [7] (same coloring). The top, middle, and lower set 
of figures correspond to Scenario I, II, and III. The "compressed wedge" corresponds to the allowed region in the upper-left of 
each plot, where the mass difference between the squark and the Higgsinos is small. The "kinematic hmit" region corresponds 
to the allowed region to the far-right of each plot, where the squark production cross section reaches the kinematic limit of the 
LHC searches. 



DISCUSSION 



We have shown that Natural Supersymmetry, vifhere 
third generation squarks decay into Higgsino-hke neu- 
trahnos and chargino, is significantly constrained by ex- 



isting LHC searches, summarized in Fig. [7j When these 
constraints are overlayed on the modifications to the 
lightest Higgs production and decay, shown in Fig. |8j 
we find distinctly different implications for the remain- 
ing allowed regions identified as the "compressed wedge" 
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{{rUq — \fi\)/mg ^ 1) and the "kinematic limit" {rriq > 
600-750 GeV). We found that the collider constraints 
arise from the totality of numerous searches at ATLAS 
and CMS that are separately sensitive, in varying de- 
grees, to squark production and decay through t — ^ tx*^, 
i — > bx^ , b — > fox", and b — )• txi ■ Our analysis incor- 
porated simulations of the signal and detector response, 
matching the experimental analyses as close as we could. 
Nevertheless, there is substantial room for improvement. 
Much of the experimental searches were designed for only 
one decay mode, or chose chargino/neutralino mass hier- 
archies that are not consistent with Natural Supersym- 
metry. We believe dedicated analyses, that take into ac- 
count the proper branching fractions and mass hierarchy, 
may well significantly improve the sensitivity. 

Natural Supersymmetry implies the wino and bino are 
sufficiently heavy that the Higgsino-like chargino and 
neutralino splittings are very small, just several GeV. 
While we focused our attention on the stop/sbottom sig- 
nals, the electroweakinos (Higgsinos) can be directly pro- 
duced at the LHC. However, the narrow splittings of the 
Higgsinos makes them extremely difficult to detect; the 
traditional search for electroweakinos is pp — )■ X2X1 ~^ 
3i + fx [139) , where the leptons come from cascade de- 
cays X2 ~^ ^XitXi ~^ Xi- As the spectrum gets 
squeezed, the intermediate W^/Z'^ go off-shell and the 
leptons they decay to are too soft to pass analysis cuts. 
Current searches are restricted to on-shell W^/Z^, so 
there is no bound from trilepton searches on degenerate 
Higgsinos. Exactly what m^o — m^o, m-± — m^o mass 

X2 Xi' Xi Xl 

splitting the experiments are sensitive to is a very inter- 
esting question, but beyond the scope of this paper. 

A more promising way to detect degenerate Higgsi- 
nos may be through monojet searches |140j . The ini- 
tial quark/anti-quark in qq — >■ xx production can emit 
hard radiation that recoils against the invisible portion 
of the event. Current monojet searches look for, among 
other physics, dark matter production after pairs of dark 
matter particles escape the detector. Re-interpreting 
the monojet bounds in terms of Higgsino pair produc- 
tion, we estimate that the current searches are sensi- 
tive to ^l ^ 100 GeV [141], competitive with the LEP 
bound on charginos |135j . With more data, this bound 
may increase, slicing into the parameter space of Natu- 
ral Supersymmetry. Additionally, there may be meth- 
ods to optimize monojet searches for Higgsinos; the ex- 
isting searches assume higher-dimensional contact op- 
erators whereas Higgsinos couple directly through elec- 
troweak gauge bosons. 

If the Natural Supersymmetry expectations for M2, 
Ml are relaxed, it is also interesting to investigate how 
sensitive the LHC will be to stop and sbottoms with light 
Higgsinos. Two effects arise from lowering M2 and/or 
Mi: the splittings between the Higgsinos increase, and 
the gaugino content of the lightest electroweakinos in- 
creases. Since stop and sbottom decays to Higgsino-like 
electroweakinos are dominated by the top Yukawa cou- 
pling, we don't anticipate significant effects on the de- 



cay branching fractions even when M2^i drop below ruq, 
opening up decays to gauginos. The larger effect is the 
increase in the splittings between the Higgsinos them- 
selves. Clearly another interesting question is to probe 
how large the splitting needs to be before the search 
strategies described here become diluted by the addi- 
tional energy from transitions between Higgsinos. Mix- 
ing the light electroweakinos with some bino, wino, or 
singlino is highly relevant for the possibility that the 
lightest neutralino could be dark matter, but this is be- 
yond the scope of this paper. 

Natural Supersymmetry may also lead to an unusual 
signal for the first and second generation squarks. One 
decay possibility for a first /second generation squark in 
Natural Supersymmetry is to a quark plus a Higgsino. 
As the first and second generation Yukawas are so small, 
the decay proceeds through the wino/bino fraction of 
the lightest neutralino and is therefore suppressed by 
0{gv'^/Mi^2)- A second decay possibility is the three- 
body decay, q j + tit or j + bib via an off-shell 
gluino. This option is suppressed by the gluino mass 
and three-body phase space, but comes with QCD cou- 
pling strength. Depending on the hierarchy of M3,M2 
as well as the mass of the light squarks relative to the 
stops/sbottoms, the three-body decay fraction can be 
substantial!^ First /second generations squark decays to 
j ti t or j bib would have several consequences that would 
be interesting to explore in more detail. Two obvious 
consequences are that the energy per final state parti- 
cle would be lower because the squarks decay to multi- 
ple particles, and the decays would contain heavy- flavor 
jets not usually associated with first/second generation 
searches. 

Finally, while Natural Supersymmetry in a low energy 
effective theory is straightforward to define and quantify, 
issues of naturalness become muddled as this is embed- 
ded into an ultraviolet completion. The obvious issue is 
the that leading-log corrections to the electroweak sym- 
metry breaking scale can quickly become a poor approx- 
imation if the renormalization group evolution is sub- 
stantial. For instance, "radiative" electroweak symmetry 
breaking arises when rn?fj^ is "driven" negative by its in- 
teraction with the stops, which clearly requires renormal- 
ization group improvement to determine the size of the 
contribution to electroweak symmetry breaking. This is 
precisely why we considered rriq and /i to be free parame- 
ters, since their separation may be much smaller than the 
leading-log approximation suggests. How this impacts 
the larger spectrum, particularly the gluino, becomes a 
highly model-dependent question. Nevertheless, we be- 
lieve our analysis has captured the essential physics of 
Natural Supersymmetry, and we remain optimistic that 



The strength of the three-body mode also depends on the mass 
character of the gluino. For Dirac gluinos the suppression in the 
three-body mode is rriq/M^ rather than I/M3, making it much 
smaller. 
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it can be discovered with continued analyses at LHC. 
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700 GeV, Ml = M2 = M3 = 2 TeV, A = 0.7, k = 0.67, 
Ax = -60 GeV, = -200 GeV and Heit = 200 GeV. 
Also, for this parameter point the low energy precision 
and flavor observables are within 2a of their experimen- 
tal valuesj^ For this point the neutralino masses are 
m^o = (197GeV, 227GeV, 416GeV, 1.98TeV, 1.99TeV) 
and the chargino masses are m^± = (200GeV, 1.98TeV). 
We checked that the branching ratios of the squarks into 
Higgsinos are within 1-2% of the an MSSM model with 
similar sfermion and Higgsino masses. 

Appendix B: Search details 



Appendix A: Realizing the observed Higgs mass and 
branching ratios of Scenario I, II and III 

In this Appendix, we consider the possibility of re- 
alizing Natural Supersymmetry scenarios with mh ^ 
125 GeV. In particular we consider the NMSSM 
model |142| where the Higgs sector of the MSSM is ex- 
tended by including a gauge singlet. The superpotential 
has the form 



Yuk 



XHuHdS + —S^ 



(Al) 



where WVuk are the usual Yukawa interactions and the 
hatted fields denote the chiral superfields. The corre- 
sponding soft supersymmetry breaking terms are 



soft 



'if. 



XA\SH,,Hd 



kA 



^S'. (A2) 



In addition to the D-term contributions to the Higgs 
mass, the additional Xv^ sin^ 2(3 contribution can help 
raise the Higgs mass above the Z-boson mass. On the 
other hand, solving the minimization conditions leads to 
the electroweak symmetry breaking condition of Eq. ([T]) 
with ^ — Aa;[^ To maximize the tree- level contributions 
to the Higgs mass, we need both the NMSSM quartic con- 
tribution as well as the usual D-term contribution, and 
thus small tan /?. Small tan /3 typically enhances the hi- 
erarchy between mn^ , "m-Ha and the electroweak scal^^ 
Due to this tension between the Higgs mass and the hi- 
erarchy of scales, we consider tan/3 S (1.5,2). However, 
we can still simultaneously realize the Natural Super- 
symmetry spectra in this paper and the observed Stan- 
dard Model Higgs mass in the NMSSM. Using NMSSM- 
tools3.2.4 |146j . we find nih — 125 GeV for the parame- 
ter space point tan/3 = 1.5, At — Af, = A^. = 0, TOj = 



For completeness, in the following we detail the im- 
portant search criteria for each collaboration's particular 
search strategy that was used in this paper. 

CMS stops, semi- leptonic, 9.7 fb"^ [43l 
Object Id: 

• jets, pt > 30 GeV, \t]j\ < 2.5, anti-kr, R = 0.5. 
Flavor tagging applied to all jets within {rijl < 2.5 

• electrons (muons), pr > 30 GeV, \t]i\ < 1.44(2.1) 

• leptons within AR = 0.4 of a jet are removed 
Basic cuts: 

• I^T > 50 GeV 

• exactly 1 lepton passing the criteria above 

• 3 or more jets, with at least one 5-tagged 
Analysis: 

• Events passing the basic selection cuts are binned 
according to the transverse mass of the MET + 
lepton system and the missing energy. Transverse 
mass is defined as 

iTT-l, jj^^^i = ^ {$tPtj - i>T ■ Pt.i) (B1) 

• The bins are, in the format {rriT^min, ^T,min)- 
(150 GeV, 100 GeV), (120 GeV, 150 GeV),' 
(120 GeV, 200 GeV), (120 GeV, 250 GeV), 

(120 GeV, 300 GeV), (120 GeV, 350 GeV),(120 GeV, 400 GeV) 

• The bins are not exclusive, so the bin with the best 
limit at a given {mt_^ , m^o ) point is used. 



The additional minimization condition of the singlet leads to a 
modified fine-tuning condition for the NMSSM. A detailed dis- 
cussion of the fine-tuning in the generalized NMSSM-like models 
can be found in Refs. l36l IT331 [Tiil . 

For alternative NMSSM scenarios utilizing large tan f}, see 
Ref. [145) 



As A and k are both somewhat large, these couplings may develop 
a landau pole before the GUT scale. The UV completion of such 
models can be realized in "fat Higgs" -like scenarios [147| . however 
a detailed study of this issue is beyond the scope of this paper. 
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ATLAS stops, semi-leptonic, 20.7 fb"^ gl] 
Object Id: 

• jets, pt > 20GeV, Irj^l < 2.5, anti-ky, R = 0.4. 
Flavor tagging applied to all jets within \ < 2.5 

• electrons (muons), px > lOGeV, \rii\ < 2.7(2.4) 

• any jet within AR = 0.2 of an electron is removed 

• subsequently, any leptons within AR ~ 0.4 of a jet 
are removed 

Basic cuts: 

• exactly 1 lepton, which must have pT^e > 25 GeV 

• 4 or more jets, at least one of which is b- 
tagged. The four hardest jets must satisfy pt > 
80 GeV, 60 GeV, 40 GeV, 25 GeV respectively. 

Analysis: 

• Most channels require top reconstruction, done as 
follows: the closest pair of jets (in AR) that satisfy 
rrijj > 60 GeV are dubbed a 'W-candidate'. This 
candidate is combined with the nearest jet (again 
in AR). For the resulting three-jet system to be 
considered a successful top-candidate, 130 GeV < 
rrijjj < 205 GeV is required. 

• Other analysis cuts include: the transverse mass 
of the — ^ system, the missing energy, the fJr- 
significance - defined as fx / \fiiT\3^i, and the 
between the missing energy (transverse) vector and 
the leading two jets. 

• In the channels designed to be sensitive to the high- 
est stop masses, two other variables are included 
amT,2 and m^jj both of which are slight variants 
on the 771^,2 variable |148j . In mT,2, as in these vari- 
ations, the visible part of the event is divided into 
two, and all partitions of the missing energy are 
scanned over. The difference between mT,2,awT,2, 
and mTj. 2 lie in whether all the visible particles are 
used, or only some of them. In amx.i, only the 
leading light jet, lepton, and highest weight 6-jet 
are taken as the visible part of the event, while in 

oiily the leading lepton and leading light jet 
are used. 

• The channels dedicated to < ^ i + are: 

1. 1 top candidate, alone with A<^-^ _^ > 0.8, 

^^i^T-j2 > f'^ > 100 GeV, fx signif. 
> 13, Mt > 60 GeV. Events passing this se- 
lection are then separated into 12 finer bins 
according to their Mt and fx- 
Mt e {60 - 90 GeV, 90 - 120 GeV, 
120- 140 GeV, > 140 GeV} 
e {100 - 125 GeV, 125 - 150, GeV 
> 150 GeV}. 



2. 1 top candidate, along with Acj)^ _. > 0.8, 

> 200 GeV, $T signif. > 13, Mt > 
140 GeV, amT,2 > 170 GeV 

3. 1 top candidate, along with > 

0.8, A0^^_^.^ > 0.8, ^T > 275 GeV, ^t 

signif. > 11, Mt > 200 GeV, amT,2 > 
175 GeV, mj, 2 > 80 GeV 

• The analysis contains three channels aimed at the 
t —i' b + final state. In these channels no top 
candidate is required. Instead there are stronger 
requirements on the pt and multiplicity of the b- 
jets, and an additional cut on rrieff, defined as the 
scalar sum of the pt of all jets with pT > 30 GeV 
summed together with the ^t magnitude and pT,i 

ATLAS stops, fully hadronic, 20.5 fb"^ ^ 
Object Id: 

• jets, Pt > 20 GeV, iTy^l < 4.5, anti-kT, R = 0.4. 
Flavor tagging applied to all jets within \T]j \ < 2.5 

• electrons (muons), pt > 10 GeV, \T]i\ < 2.7(2.4) 

• any jet within AR = 0.2 of an electron is removed 

• subsequently, any leptons within AR — 0.4 of a jet 
are removed 

Basic cuts: 

• zero leptons passing the above criteria 

• ^T> 130 GeV 

• 6 or more jets, where jets satisfy pT > 33 GeV, 
\r]j\ < 2.8. The leading two jets must have pt > 
80 GeV, and at least two jets are 6-tagged. 

Analysis: 

• 2 three-jet clusters are formed from from the list 
of jets as follows: the three jets that are closest 
in the (f> — i] plane are taken as one such cluster, 
removed from the list, then the process is repeated 
to extract the second group. The mass of these 
three-jet clusters is required to lie within 80 GeV < 
rrijjj < 270 GeV in order to select events with two 
hadronic tops in the final state. 

• The transverse mass of the ]pT~b system, where the 
b closest in Acj) is used is required to be > 175 GeV 
to remove leptonic tt background 

• _. > 0.27r, where A(j)^ . is the angle between 

the missing energy vector and the closest jet. This 
cut is designed to remove backgrounds from mis- 
measured jets. 
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• The remaining events are binned according to ]pT'- 
$T > 200 GeV, > 300 GeV and > 250 GeV. 
Only the strongest limit at a given (to^^ , m^o ) point 
is used. 



CMS sbottoms, multi-& + ^t, 11.7 fb"^ [35] 
Object Id: 

• jets, pT > 50Ge\pJ |r;j| < 3.0, anti-ky, R = 0.5. 
Flavor tagging applied to all jets within \7]j \ < 2.5 

• electrons (niuons), pr > 10 GeV, \T]e\ < 2.4(2.1) 

• any jet within AR = 0.4 of a lepton is removed 
Basic cuts: 

• zero leptons 

• at least 2 jets. The hardest jet must lie within 
|?7| < 2.5 and the leading two jets must have px > 
twice the nominal jet px requirement. Nominally 
this is > 100 GeV but for events with \ow-Ht this 
cut may be softer. Events with high-p-r jets (i.e. 
passing nominal jet criteria) at \ri\ > 3.0 are vetoed. 

• Ht > 275 GeV, where Ht is the scalar sum of the 
Pt of all jets in the events. 

Analysis cuts: 

• All visible objects in the event are grouped into 
two mega-jets, following the criteria given in Ref. []. 
The degree to which the two megajets balance each 
other is captured by the variable ar = Mt ^ - ' 
fraction of the transverse energy of the subleading 
(in Pt) megajet relative to the transverse mass of 
the megajet pair. Requiring ut > 0.55 greatly sup- 
presses multijet QCD backgrounds. 

• Events surviving the ax cut are categorized accord- 
ing to the jet and 6-jet multiplicities, then binned 
in Ht- 

• The jet multiplicity categories are Njet =2 — 3, 
and Njet = 4"*". Within each jet multiplicity cat- 
egory, Nb = 0,1,2,3,(4) is considered (obviously 
Nb = 4: is only considered in the Njet >= 4 
class). For a given {Njet,Nb), the Ht is binned as 
[275 - 325 GeV], [325 - 375 GeV], [375 - 475 GeV]„ 
[475 - 575 GeV], [575 - 675 GeV], [675 - 775 GeV], 
[775 - 875 GeV] and > 875 Ge\f5 



This requirement is scaled down to 37 GeV, or 43 GeV for events 
with \ow-Ht- 

For the samples with the highest 6-multiplicity, only three hx 
bins are used, [275 - 325 GeV], [325 - 375 GeV], > 375 GeV. 



• For the direct sbottom search, only the Njet = 2 — 
3, iVf, = 0, 1 categories are used to set limits. Since 
the Ht bins are orthogonal, all Ht bins across both 
categories are taken together to form a combined 
limit. 

ATLAS sbottoms, multi-6 + fJr, 12.8fb"^ [44] 

Object Id: 

• jets, PT > 20 GeV, \rij\ < 2.8, anti-kr, R = 0.4. 
Flavor tagging applied to all jets within \rij \ < 2.5 

• electrons (muons), pT > 10 GeV, jry^l < 2.7(2.4) 

• any jet within AR ~ 0.2 of an electron is removed 

• subsequently, any leptons within AR — 0.4 of a jet 
are removed 

Basic cuts: 

• zero leptons 

• 2 or more jets, with 2 or more &-tags 

• > 150 GeV 

Analysis: 

• After basic selection, 3 event categories are set 
up, each with slightly different requirements. The 
categories are not exclusive: 

1. leading jet pt > 150 GeV, subleading jet 
Pt > 50 GeV, no other jets with pT > 50 GeV. 
Both the leading two jets must be tagged as b 
jets. To reduce multijet QCD, A(j)^^_^^ > 0.4 

and ^T/iTieff > 0.25 are required. Here 
rueff is the scalar sum of the missing en- 
ergy and the pt of the hardest three jets sat- 
isfying basic jet requirements (meaning they 
must be harder than 20 GeV only) ,mef / = 

^T + J2'i=i PT,ji ■ Within this category, events 
are further binned according to their contra- 
transverse mass, see Ref. |149| for definition. 

2. similar to category 1.) but the pt require- 
ments are adjusted to > 200 GeV, > 60 GeV 
for the leading and subleading jets. The lead- 
ing two jets still must be flavor tagged, and the 
A(j)^ and the ^T/rUeff are unchanged. 

3. More than 2 jets are required with the lead- 
ing jet having pt > 130 GeV. The two sub- 
leading jet must have pt > 30 GeV, but be- 
low 110 GeV. Unlike the previous categories, 
the first two categories, the leading jet is not 
required to be a 6-jet. Instead the leading 
jet must have light flavor (it is anti-tagged). 
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while the subleading two jets must be tagged 
as &-jets. The /S.(j)^^_.^ and the JfiT/meff 

are unchanged, but there is an additional re- 
quirement that the scalar sum of the pr of 



all jets beyond the leading three is small, 
< 50GeV. This category is divided into two 
subcategories with different and PT,ji re- 
quirements. 
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